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ABSTRACT 

We present a survey of mid-infrared gas-phase lines toward a sample of 76 circumstellar disks around low 
mass pre-main sequence stars from the Spitzer "Cores to Disks" legacy program. We report the first detections 
of [Ne II] and [Fe I] toward classical T Tauri stars in ~ 20 % respectively ~ 9 % of our sources. The observed 
[Ne II] line fluxes and upper limits are consistent with [Ne II] excitation in an X-ray irradiated disk around stars 



with X-ray luminosities L x = 10 29 - 10 31 ergs _1 . [Fel] is detected at ~ 10" 5 - 10" 4 L Q , but no [SI] or [Fell] 
is detected down to ~ 10~ 6 L©. The [Fel] detections indicate the presence of gas-rich disks with masses of 
> O.lMj. No H2 0-0S(0) and S(l) disk emission is detected, except for S(l) toward one source. These data 
give upper limits on the warm (T ~ 100- 200 K) gas mass of a few Jovian masses, consistent with recent 
T Tauri disk models which include gas heating by stellar radiation. Compact disk emission of hot (T > 500 K) 
gas is observed through the H2 0-0S(2) and/or S(3) lines toward ~ 8 % of our sources. The line fluxes are, 
however, higher by more than an order of magnitude than those predicted by recent disk models, even when 
X-ray and excess UV radiation are included. Similarly the [Ne II]/H2 0-0 S(2) ratios for these sources are lower 
than predicted, consistent with the presence of an additional hot molecular gas component not included in 
current disk models. Oblique shocks of stellar winds interacting with the disk can explain many aspects of the 
hot gas emission, but are inconsistent with the non-detection of [S I] and [Fe II] lines. 

Subject headings: infrared: ISM - planetary systems: protoplanetary disks 
stars: evolution - stars: low-mass 



survey - circumstellar matter - 



1. INTRODUCTION 

Circumstellar disks are a natural and important phe- 
nomenon in low-mass star formation. The ability to trace the 
evolution of the disk dust and gas content is crucial to un- 
derstanding their chemistry and physics and the formation of 
planets. Observed disks reveal a large source to source vari- 
ation and suggest a complex evolution from young gas-rich 
disks to tenuous debris disks. Disk geometries are observed 
to range from almost flat to strongly flaring disks (Dulle- 
mond & Dominik 2004), and disks with large inner cavities 
are found (Calvet et al. 2002; Bouwman et al. 2003; Forrest 
et al. 2004; Brown et al. in prep.). Observations of silicates 
and spectral energy distributions (SEDs) present evidence for 
grain growth and settling of large grains to the disk midplane 
(van Boekel et al. 2003; D' Alessio et al. 2006; Kessler-Silacci 
et al. 2006). Different degrees of grain heating and radial mix- 
ing in the disks (van Boekel et al. 2005), and varying PAH 
abundances affect the heating of the upper layers of the disk 
(Habart et al. 2006; Geers et al. 2006). Most of these disk 
properties are derived from infrared solid-state features and 
from near-infrared to millimeter continuum observations and 
associated SED modeling. However, such data give little in- 
formation about the gas in the disk. 

Gas plays an important role in the structure and evolution of 
disks, including the temperature and density gradients, chem- 
istry, dust dynamics, and eventually the formation of gas-rich 
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and gas-poor planets (Gorti & Hollenbach 2004). Observa- 
tional diagnostics of the physical conditions of the gas, in 
particular its mass and temperature, are therefore highly rel- 
evant to studies of disk evolution and planet formation. Cen- 
tral questions are how the gas in the disk is dissipated, what 
drives the disk heating and gas clearing, and what the asso- 
ciated timescales are. The small number of observed transi- 
tional objects between the classical T Tauri phase (large Ha 
equivalent width, thought to be accreting) or the weak-line 
T Tauri phase (small Ha, non accreting), both with massive 
optically thick disks, and the more evolved phase with opti- 
cally thin or no disks, suggests that disk clearing timescales 
are short (fewx 10 5 yr) compared to the disk lifetime of a few 
Myr (Cieza et al in prep.; Haisch et al. 2001). In addition, 
combined near-infrared (IR), mid-IR, and millimeter obser- 
vations imply that the disk clearing happens nearly simulta- 
neously across the disk (see e.g. Kenyon & Hartmann 1995; 
Hartmann et al. 2005; Takeuchi et al. 2005). Alexander et al. 
(2006) present a new evolutionary model combining viscous 
evolution with photoevaporation of the disk to address these 
issues. In this mechanism, the disk is cleared through an evap- 
orative flow originating from the disk surface layers as a re- 
sult of UV and/or X-ray heating (see e.g. Hollenbach et al. 
2000; Clarke et al. 2001; Kamp & Dullemond 2004; Dulle- 
mond et al. 2007). Information about the temperature and 
mass of the warm gas and constraints on the details of the 
heating processes across the disk are therefore vital for ad- 
vancing the current models. 

Models of gas heating include UV and X-ray radiation from 
the star itself as well as possible excess radiation due to ac- 
cretion (e.g. Jonkheid et al. 2004; Gorti & Hollenbach 2004; 
Kamp & Dullemond 2004; Nomura & Millar 2005; Jonkheid 
et al. 2007). The resulting gas temperatures in the surface 
layers out to large radii are significantly higher than those of 
the dust as a result of the photoelectric effect on small grains 
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and PAHs. Gas temperatures may reach values of up to a few 
thousand K. Once grains have grown to /im sizes or larger, 
however, the gas temperature is significantly decreased unless 
PAHs are still present (Jonkheid et al. 2004, 2006, 2007). No- 
mura & Millar (2005) include explicitly the excitation of H2 
by UV and collisions for a disk around a T Tauri star with 
and without excess UV and make predictions for line intensi- 
ties which can be tested against observations. Glassgold et al. 
(2007) have studied the case of X-ray heating and have shown 
that the ionized neon fine-structure line emission may pro- 
vide unique tracers of X-ray heating in the disk surface since 
neon cannot be photoionized by radiation with energies below 
13.6 eV. 

Observational studies of the gas and its temperature have 
mostly focused on the hot inner and the cold outer regions 
of disks. High-resolution CO v=l-0 vibration-rotation lines 
at 4.7 fim (Najita et al. 2003; Brittain et al. 2003; Blake & 
Boogert 2004) and H 2 2. 1 fim 1-0 S(l) (Bary et al. 2003) show 
gas with temperatures ~ 1000- 3000 K in the surface layers 
out to - 1 AU. H 2 emission from SVS 13 (Carr et al. 2004) 
and molecular absorption of C2H2, HCN, and CO2 in the disk 
of IRS 46 (Lahuis et al. 2006a) also indicate hot temperatures 
in the inner few AU of several hundred K. In contrast, mil- 
limeter CO surveys probe the cold gas throughout the outer 
disk where the dust is optically thin (e.g. Koerner & Sargent 
1995; Duvert et al. 2000; Thi et al. 2001; Dutrey et al. 2003; 
Dent et al. 2005). Its use as a gas mass tracer is however 
limited as a result of both strong photodissociation at low ex- 
tinction and freeze-out in the disk interior. Gas temperature 
determinations range from <20 K near the midplane (Dartois 
et al. 2003; Pietu et al. 2007) to 40 K or higher in the interme- 
diate and surface layers (van Zadelhoff et al. 2001; Qi et al. 
2006). 

The mid-IR H2 and atomic fine structure lines are best 
suited as direct tracers of the warm (^100 K) gas in the in- 
termediate zones of disks at radii of a few AU out to sev- 
eral tens of AU, i.e., the planet-forming zones of disks. The 
Infrared Space Observatory (ISO) provided the first opportu- 
nity to probe this warm gas in disks around Herbig Ae and 
T Tauri stars. Thi et al. (2001) suggested that large amounts 
(<~ 0.01 - 100My) of gas could reside in disks around young 
T Tauri stars, but this has not been confirmed by subsequent 
ground-based observations (Richter et al. 2002; Sheret et al. 
2003; Sako et al. 2005). 

The sensitive InfraRed Spectrograph (IRS) (Houck et al. 
2004) on board the Spitzer Space Telescope (Werner et al. 
2004) brings the detection of these lines within reach for 
young solar mass stars in nearby star forming regions. The 
combination of high sensitivity, moderate spectral resolution 
R = A/ A A = 600, and modest spatial resolution makes Spitzer 
well suited for the direct study of the gas in and around low- 
mass young stars in nearby (< 300 pc) clouds through the 
mid-IR lines of various species. 

We present here an overview of gas-phase lines detected in 
disks observed in the Spitzer legacy program "From Molecu- 
lar Cores to Planet Forming Disks" ("Cores to Disks" or c2d) 
(Evans et al. 2003), which has collected a large sample of 
IRS spectra toward sources in the nearby Chamaeleon, Lupus, 
Perseus, Ophiuchus, and Serpens star-forming regions. High- 
S/N 5-38 [im spectra have been obtained for 226 sources at 
all phases of star and planet formation up to ages of ^5 Myr. 
From this sample, 76 disk sources, identified by showing ei- 
ther the 10 or 20 fim silicate bands in emission, have been 
selected. In Sections 2 and 3 the source selection and data 



reduction are explained. In Section 4 the observed atomic 
fine-structure and H2 emission lines and the derived parame- 
ters are presented. In Section 5 the results are reviewed in the 
context of currently available disk models. This paper forms 
a complement to the searches for the mid-infrared lines of H2 
and other species toward more evolved disks studied in other 
Spitzer programs (e.g. Hollenbach et al. 2005; Pascucci et al. 
2006; Pascucci et al. 2007). 

2. OBSERVATIONS 

The data presented in this study were selected from the 
sample of IRS spectra observed within the Spitzer c2d legacy 
program. The c2d IRS program consists of two programs of 
comparable size, referred to as the first- and second-look pro- 
grams. The first-look program (PID#172) was restricted pri- 
marily to known low-mass young stars, embedded YSOs and 
pre-main-sequence stars with disks with masses M < 2 Mq 
and ages < 5 Myr, and a sample of background stars. A few 
Herbig Ae stars are included as well. The c2d source se- 
lection criteria were defined to be complementary to those 
of the Spitzer legacy program "The Formation and Evolu- 
tion of Planetary Systems" (FEPS, Meyer et al. 2002). The 
second-look program (PID#179) was, for the most part, de- 
voted to IRS follow-up spectroscopy of sources discovered 
in the IRAC and MIPS mapping surveys, including a newly 
discovered cluster of young stars in Serpens (Harvey et al. 
2006). For all first-look observations, the integration times 
for the short-high (SH) and long-high (LH) modules (R = 600, 
10-37 fim) were chosen such that theoretical signal to noise 
ratios (S/Ns) of at least 100 and 50 on the continuum were 
obtained for sources brighter and fainter than 500 mJy, re- 
spectively. Deeper integrations were not feasible within the 
c2d program. Spectra taken using the short-low (SL) or long- 
low (LL) modules (R = 60- 120, 5-14 /xm and 14-38 fim re- 
spectively) always reach theoretical S/N ratios greater than 
100. For the second-look IRS targets similar S/N limits were 
obtained wherever possible. However, since the second-look 
contained a number of very weak sources (down to a few mJy) 
this was not always achieved. 

2.1. Source selection 

The sources presented in this paper were all selected to 
show either of the 10 /jm or 20 /im silicate bands in emission. 
A total of 76 sources were chosen; see Kessler-Silacci et al. 
(2006) for the 47 first-look disk sources with silicate emission. 
This selection excludes most edge-on disk sources (;' > 65 de- 
grees) such as CRBR 2422.8-3423 (Pontoppidan et al. 2005), 
IRS 46 (Lahuis et al. 2006a) and the 'Flying Saucer' (Pon- 
toppidan et al. 2007), with the exception of the high incli- 
nation sources EC82 (Pontoppidan & Dullemond 2005) and 
VV Ser which are included in this paper. Gas-phase lines to- 
ward edge-on disk sources will be discussed in a separate pa- 
per together with the embedded class and I sources (Lahuis 
et al, in prep.). The selected sources are listed in Table 1 
which gives the basic observing and source parameters, e.g. 
the adopted distances. 

2.2. SH mini maps 

In an early phase of the c2d project molecular hydro- 
gen and [Nell] lines were detected toward some of the c2d 
sources. As part of the second-look program, follow-up mini- 
maps were taken using the SH module to check for extended 
emission at offsets positions of ^10-15" with respect to the 
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Fig. 1 . — Observations of Hi 0-0 S(l), S(2), and [Nell] emission observed on and off source toward Sz 102, EC 74, and Haro 1-17. The on 
source observations were observed in the c2d first-look program, the follow up off source observations in the second look program (see Sec. 2). 
The images show Spitzer composites of IRAC1 (blue), IRAC2 (green), and IRAC4 (red) with the sources in the middle. The gray lines in the 
spectra show the total (compact source plus extended component) observed emission. The black lines the source emission after correction for 
the extended component. The red lines show Gaussian fits to the compact source emission if observed at 3a or more. Although the S/N in the 
off source positions is often poor, it illustrates the problem of extended and non-uniform H2 emission. The [Ne II] is almost always limited to 
the source itself. 
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Fig. 2. — Illustration of the Spitzer IRS cross dispersion profile used in the optimal extraction (see Sec. 3.1). The left plot shows a fit to the IRS 
SH order 1 1 data of GW Lup, a source with a moderate but clear sky component in the IRS spectra. The rsc data (bed data before flatfielding) 
of both dither position (black plusses) is shown, normalized, collapsed along the dispersion direction, and corrected for the cross dispersion 
dither offsets. Overplotted is the combined fit of the source profile plus the extended emission in gray and the extended emission in dark gray. 
The shape of the extended emission reflects the IRS Hatfield of the (for this source) assumed uniform extended emission. The right plot shows 
a comparison of an IRS PSF profile (Sine plus harmonics) compared to the profiles of an undistorted Sine, a Gaussian, and a Lorentzian profile 
with the same FWHM. Note the significant variation in the strength and shape of the profile wings. The correct characterization of both the 
width and the wings of the profile for all IRS orders is essential for extracting the proper source and sky spectra. 



sources. Five maps were denned to include off source obser- 
vations around eight sources. Three of these are disk sources 
included in our sample, Sz 102 (Krautter's star), Haro 1-17, 
and EC74. Figure 1 shows the observed H 2 0-0 S(l), S(2), and 
[Ne II] emission from the first-look on source observations to- 
gether with the off source observations from the second-look 
mini maps. The maps show that most of the H2 emission is 
extended, especially for the S(l) line. However, fine-structure 
lines are usually seen to be limited to the source itself. Subse- 
quent results using the c2d optimal extraction procedure (see 
Section 3.1) confirm the conclusions drawn from the analysis 
of the mini-maps. 

The SH maps are not complete since the time allocated 
for c2d spectroscopy did not allow to observe fully sampled 
maps. The prime purpose of the maps is to confirm the pres- 
ence or absence of extended emission. More extended map- 
ping observations will be required to study the detailed spatial 
distribution and extent of the large scale emission component. 

3. DATA REDUCTION 

The c2d reduction pipeline (Lahuis et al. 2006b) was used 
to reduce the IRS data, starting from the S13 and S14 archive 
data. The same c2d pipeline products as included in the final 
c2d Legacy data delivery 5 were used for the spectral line anal- 
ysis (see Sec. 3.3). Most of the analysis focused on the SH 
and LH data, since the SL and LL data are generally limited 
by the low line/continuum ratio. The SL data were included 
and used to search for higher-excitation H2 0-0 lines, in par- 
ticular the S(3) transition. 

3.1. Separating disk and cloud emission - optimal 
extraction 

A major concern when studying emission lines from young 
stellar objects (YSOs) is the possible contribution of extended 
(envelope or local cloud) emission in the sometimes complex 
star forming regions. The spatial distribution of the emission, 

5 The c2d legacy data are accessible at 
http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html. 



both in the continuum and in spectral lines, often prohibits the 
use of 'sky' observations alone to correct for extended emis- 
sion components. For this reason the c2d team has developed 
an optimal extraction algorithm for IRS pointed observations. 

The Spitzer diffraction limited beam is ^4-5 arcsec for 
the SH module (10-19.5 fjm) and - 7- 10 arcsec for the LH 
module (19-37 /im). At a distance of 100 parsec this corre- 
sponds to physical sizes of -400-500 AU and ~700-1000 
AU, respectively. The clouds observed in the c2d program are 
located at distances ranging from 125 parsec (Ophiuchus) to 
260 parsec (Serpens) increasing the physical area observed. 
The full IRS aperture in the SH and LH spatial direction is 
^2.5-3 times larger than the beam size. At the observed 
cloud distances this means that the aperture probes physical 
scales of several thousand AU. This makes it ideally suited 
for detecting cold or shocked H2 emission from the extended 
(remnant) envelope, outflows, or the diffuse local cloud emis- 
sion. For the disk sources studied in this work however the lo- 
cal cloud emission will potentially confuse the compact disk 
emission. Distinguishing between compact (disk) and ex- 
tended (remnant envelope, outflow, or diffuse cloud) emis- 
sion is therefore of vital importance for studying the emission 
lines originating in the circumstellar disks. The optimal PSF 
extraction developed by the Spitzer c2d legacy team (Lahuis 
et al. 2006b) allows separation of the two components for all 
sources. The mini-maps (see Section 2.2) observed around 
selected sources confirm the results of the optimal extraction. 
See Figure 1 for an example of extended H2 emission but 
compact [Nell] emission. Other examples include Figures 
3 and 4 in Geers et al. (2006) for separating extended cloud 
and compact disk PAH emission. 

The optimal extraction uses an analytical cross-dispersion 
point spread function (PSF) for the source profile plus an ex- 
tended emission component, to fit the observed crosstalk or 
straylight corrected echelle (SH and LH) and longslit (SL and 
LL) images. The PSF is described by a sine function with 
a harmonic distortion component which results in a signifi- 
cant broadening of the profile wings (see Fig. 2). The wave- 
length dependence of the PSF parameters, the order trace, the 
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width, and the harmonic distortion, are characterized using a 
suite of high S/N calibrator stars. For the extended emission 
component the flatfield cross-dispersion profile is used. The 
flux calibration derived from the calibrator stars using Cohen 
templates and MARCS models (Decin et al. 2004) provided 
through the Spitzer Science Center. Lahuis et al. (2006b) give 
more details about the characterization and calibration of the 
c2d optimal extraction. The optimal extraction returns the to- 
tal flux (the source flux plus extended emission in the IRS 
beam) plus an estimate of the extended emission component. 
Error estimates are derived for both the total emission and 
the extended emission component. The S/N of the extended 
emission component can vary significantly depending on the 
quality of the raw image data and on deviations of the ex- 
tended emission from the assumed uniformity across the IRS 
slit. Therefore, care has to be taken when subtracting the ex- 
tended emission from the total flux signal to retrieve the com- 
pact source emission. In some cases, a fit to the extended 
continuum and line emission is used to avoid adding in sur- 
plus noise from the extended emission component. The un- 
certainty on the fit to the extended emission is propagated into 
the error of the compact source signal. 

3.2. 1-D spectra 

After extraction, the 1-D spectra are corrected for in- 
strumental fringe residuals (Lahuis & Boogert 2003), order 
matching is applied, and a pointing flux-loss correction is per- 
formed to the compact source component. Pointing offsets 
up to a few arcsec can have a noticeable impact on the de- 
rived fluxes of lines observed with the SH and SL modules, 
e.g. H 2 0-0S(1), S(2), [Nell], and [Nelll]. For example, dis- 
persion offsets within the nominal 3a pointing uncertainty of 
Spitzer (~ 1" for medium accuracy peakup) can lead to SL 
and SH flux losses up to ~ 10% depending on wavelength. 
For all targets, a combination of either the SH, LH, and SL, 
SH and LH, or SH and SL modules is available. This allows 
correction of the pointing related flux losses with an accuracy 
given by the S/N of the data in the module overlap areas. A 
detailed description of the c2d pipeline (including extraction, 
defringing, pointing flux loss correction) and the c2d legacy 
products is given in the "c2d Spectroscopy Explanatory Sup- 
plement" (Lahuis et al. 2006b). 

3.3. Spectral Analysis 

The SH and LH modules of the IRS instrument cover the 
positions of the three lowest pure rotational lines of H2 and 
emission of [Nell] (12.8 /an), [Nelll] (15.55 /zm), [Fel] (24 
Aim), [Fell] (17.9 and 26.0 urn), [SI] (25.25 urn), [Sill] (18.7 
/xm), and [Si II](34.8 /an) at a resolving power of R= A/AA = 
600. Line fitting and flux integration is done using routines 
from OSIA 6 . 

As discussed in Section 3.1, the extended emission compo- 
nent, both in the continuum and the spectral line, is subtracted 
from the spectrum prior to line fitting. Uncertainty estimates, 
as listed in Table 2, are derived from the residuals after line fit- 
ting, or, in the absence of a spectral line, using the line width 
derived from the instrumental resolution. The uncertainty de- 
rived from the extended emission is added into the uncertainty 
estimate of the source component. As a result, the 1-a uncer- 
tainty estimates can vary widely for sources with a similar 

6 OSIA is a joint development of the ISO-SWS consortium. Contribut- 
ing institutes are SRON, MPE, KUL and the ESA Astrophysics Division. 
http://sws.ster.kuleuven.ac.be/osia/ 
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Fig. 3. — Detections of neon lines at the 3a level or more toward 
the c2d sample of T Tauri stars with disks. Of all the H2 and atomic 
fine structure lines [Nell] is most convincingly detected toward ~ 
20% of the sources. Plotted in gray is the total observed emission 
(compact source + extended component) and in black the emission 
after correction of the estimated sky component. None of the sources 
show a significant extended component. 



continuum flux and integration time. This may for example 
be the result of the presence of artifacts resulting from hot 
pixels or small variations in the extended emission which are 
not accounted for in the spectral extraction which assumes a 
constant extended emission component. 

Typical mean 3a uncertainties prior to subtraction of the 
extended component for the high resolution modules range 
from ~ 1 x 10~ 16 -2 x 10~ 15 erg cm" 2 s" 1 with positive and 
negative extremes of ~ 5 x 10~ 17 ergcm~ 2 s _1 and ~ lx 
10~ 14 erg cm" 2 s" 1 . The uncertainties are comparable to those 
from Pascucci et al. (2006) for FEPS observations using on 
source integration times similar to those used for the c2d sam- 
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Fig. 4. — Observed line strengths (hatched bars) and upper limits 
(solid gray bars) of the major atomic lines. Included with vertical 
lines are model predicted line strengths for [Nell] and [Nelll] by 
Glassgold et al. (2007). The model incorporates UV and X-ray heat- 
ing of the circumstellar disk. The two lines represent two extreme 
thermal models, X-ray heating dominant (dotted line) or mechanical 
heating dominant (dashed line). 

pie. 

4. RESULTS 
4.1. Atomic fine-structure lines: Neon 

Of all the atomic fine structure lines covered by the SH and 
LH modules, the [Nell] 12.8/im transition is most convinc- 
ingly detected and shows the strongest source emission. The 
higher excitation [Nelll] 15.5/im line is tentatively detected 
toward Sz 102 but not toward any of the others sources in our 
sample. Figure 3 shows all [Ne II] and [Ne III] lines detected 
at 3(7 or more; see also Fig. 1 of Geers et al. (2006) for the 
[Nell] 12.8 /im line toward T Cha. Plotted in Figure 3 is the 
total continuum-subtracted observed [Ne II] emission in gray 
and the compact source emission after correction for extended 
line emission in black. Taken together, [Ne II] emission is ob- 
served in the spectra of 15 T Tauri sources (~ 20% of the 
sample). These are the first reported detections of [Nell] to- 
ward disks around classical T Tauri stars. The optimal extrac- 
tion method, together with the limited mini-maps (Sect. 2.2 
and Fig. 1), show that the emission is indeed associated with 
the source itself. 

The observed line fluxes and upper limits of [Nell] and 
[Nelll] are listed in columns 8 and 9 of Table 2. When 
line fluxes are compared (observed or with models) the line 
strength is converted to solar luminosities since the sample 
is observed toward sources from multiple clouds and com- 
pared to model predictions using different assumed distances. 
Figure 4 shows the distribution of the observed line strengths 
(hatched bars) and upper limits (solid gray bars). 

4.2. Atomic fine-structure lines: Other species 

[Fe I] at 24/im is the only other species besides [Ne II] with 
clear detections toward seven sources (~ 9 % of the sample). 
Of the other atomic lines there is one detection of [Fe II] at 
18/im and one of [S III] at 18.7 /im, in different sources. [S I] 
at 25 /im, [Fell] at 26/im, and [Sill] at 34 /im are not de- 
tected. The derived line fluxes and upper limits are listed in 
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columns 9- 14 in Table 2. The detected lines are plotted in 
Figure 5 whereas Figure 4 shows the distribution of the ob- 
served line strengths. 

4.3. Molecular Hydrogen 

H2 emission is detected toward a small number of sources. 
Figure 6 shows the observed H2 0-0S(2) and S(3) emission 
lines detected at 3a or more. The total (compact and extended 
line emission) observed H2 emission (in gray) and the emis- 
sion after subtraction of the extended line emission (in black) 
are shown. The H 2 0-0S(0) and S(l) lines are seen toward 
some sources such as HD 132947 and Sz 102 (Fig. 1) but they 
are found to be mostly extended. Only for Sz 102 compact 
source emission of H2 0-0 S(l) is tentatively detected. Neither 
S(0) nor S(l) are seen toward HD 135344 and HD 163296. 
The 3cr upper limits for HD 135344 are a factor of ~ 2- 3 
lower than the tentative detections in Thi et al. (2001), while 
for HD 163296 the upper limits are comparable to the ISO 
SWS line fluxes. The observed line fluxes and upper limits of 
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the H 2 0-0 S(0), S(l), S(2), and S(3) emission lines are listed 
in columns 4 - 7 of Table 2 . Figure 7 shows the distribution of 
observed line strengths (hatched bars) and upper limits (solid 
gray bars) of H 2 0-0S(0), S(l), S(2), and S(3). 

In the simplest analysis, the H 2 excitation is assumed to 
be in local thermal equilibrium (LTE) (e.g., Thi et al. 2001) 
with an ortho-to-para ratio determined by the kinetic temper- 
ature of the gas (following Sternberg & Neufeld 1999). For 
gas temperatures 100, 150, and >200K, the ortho-to-para ra- 
tios are 1.6, 2.5, and 3, respectively. Assuming optically thin 
emission, the integrated flux of a rotational line J u — > Ji for a 
given temperature T ex is 
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where A is the wavelength of the transition, N(H 2 ) the to- 
tal column density, A u i the spontaneous transition probabil- 
ity, and f2 the source size. For high enough densities (n > 
10 3 cm" 3 ), the population x u follows the Boltzmann law 
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Fig. 7. — Observed line strengths (hatched bars) and upper limits 
(solid gray bars). The vertical lines present model T Tauri disk line 
strengths of H 2 S(0), S(l), S(2), and S(3) from Nomura et al. (2007). 
The models incorporate UV and X-ray heating of the circumstellar 
disk for three grain size distributions (a max = lOpim, 1mm, and 10cm 
represented by dotted, dashed, and dashed-dotted lines respectively). 
For S(0) the line strength increases with increasing maximum grain 
size, while for S(2) and S(3) the line strength decreases as the maxi- 
mum grain size increases. 



where Ej is the energy of the upper level, is the nuclear 
statistical weight (1 for para and 3 for ortho H 2 ), and Q(T ex ) 
the partition function for the given excitation temperature T ex . 

Using the above equations, excitation temperatures, column 
densities and H 2 gas masses can be derived from the observed 
line fluxes and upper limits. If either S(0) or S(l) are detected 
an upper or lower limit on the temperature of the warm gas is 
derived, but if neither are detected a temperature of 100 K is 
assumed for the warm gas. If two or more higher excitation 
lines (S(2) and higher) are detected a temperature for the hot 
component is derived, while if no or only one of the higher ex- 
citation lines is detected a temperature of 1000 K is assumed. 
For Sz 102 and EC74 temperatures of r hot ~ 700-800 K could 
be found for the hot component. For all other sources no tem- 
peratures could be derived for either component. 

The column density averaged over the IRS aperture can be 
derived from the above equations, given the distance to the 
source. For all sources in our sample the emitting source 
size in the disk is smaller than the IRS aperture (Sec. 3.1) 
and since this is unknown a typical emitting disk region is 
assumed. For the warm component a source with a radius 
r = 100 AU is assumed and for the hot component a source 
with a radius r = 2 AU. The derived or assumed temperature 
plus the (upper level) column densities give a total column 
density, which in turn gives the total H 2 gas mass in Jovian 
masses, M = irr 2 xN x 2m H /Mj with m H = 1 .674 • 10" 24 gr 
andMj = 1.9-10 30 gr. 

The derived H 2 parameters for both the warm and hot gas 
component are listed in Table 3. Figure 8 shows the distribu- 
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tion of the derived H2 masses for the assumed temperatures of 
1 00 K and 1 000 K, respectively. 



4.4. Correlations 

Figure 9 shows the observed line strengths and upper limits 
of H 2 0-0S(2) [Nell], and [Fel] as functions of the mid-IR 
luminosities (in Solar luminosities), the 12.8-26/im spectral 
index, and the effective temperature. The mid-IR luminosity 
is integrated using the 12.8 and 15.5^m continuum points de- 
rived in the [Nell] (12.8/zm) and [Nelll] (15.5/um) spectral 
line fits. The effective temperature is derived from the stellar 
type (Gray & Corbally 1994; Baraffe & Chabrier 1996). In 
all panels the sources with H2 0-0S(2), [Nell], or [Fel] de- 
tections are plotted with large symbols whereas the sources 
without line detections are plotted with small gray symbols. 

The line strength plots as function of the mid-IR luminosity 
clearly show the detection limits of the Spitzer IRS instrument 
with an obvious increase of the upper limits with increasing 
continuum. Therfore no definite conclusions can be drawn 
about the apparent correlations, and the detection rate could 
be increased significantly with higher S/N for strong contin- 
uum sources. However the plots do show that there are no low 
luminosity sources with strong H2 0-0S(2) or [Nell] emis- 
sion. No correlation with the 12.8 -26^m spectral index is 
observed. The detections are distributed over a range of spec- 
tral indices, illustrating that the line detections are not limited 
to a single type of disk source (e.g., flat or flaring). 

The correlation with the effective temperature shows a dif- 
ferentiation between [Ne II] and [Fe I] compared with H 2 0- 
0S(2). The [Nell] and[FeI] line strengths show a similar 
correlation with effective temperature as with the mid-IR lu- 
minosity, but with more scatter. Also, upper limits are seen 
below the correlation line for sources with detections. There 
are a few detected sources deviating from the observed trend, 
such as the cold disk source T Cha (see Brown et al. in prep.) 
detected in [Ne II] and the Herbig Ae star HD 132947 detected 
in [Fel]. TCha is located at a distance of 66 pc, much closer 
than the majority of sources in our sample. Sources with a 
similar [Ne II] line strength as T Cha but at the distances of 
the nearest star-forming clouds would go undetected at the 
sensitivity limits of the current sample. For HD 132947 the 
distance is unknown and the assumed distance of 60 pc is the 
lower limit from Tycho. The H 2 0-0 S(2) line, although de- 
tected for only a small number of sources, differs from [Ne II] 
and [Fe I] in that all sources are concentrated around a single 
effective temperature. This may be real, but it could also be 
the result of a S/N selection. More sensitive observations will 
be required to draw firm conclusions. 

5. DISCUSSION 

The most significant outcome of this survey is the detection 
of [Nell] emission toward 15 T Tauri stars, ^20% of the 
sample. [Fe I] is also seen toward 7 sources, ~ 9 % of the 
sample (see Sect. 4.1 and 4.2). 

5.1. [Nell] 

Since neon cannot be ionized by photons with energies less 
than 21.4 eV (and Ne + 41.0eV), the detection of [Nell] is 
evidence for the presence of higher energy photons in the cir- 
cumstellar environment, specifically EUV photons or X-rays 
originating from either the stellar chromosphere or (accretion) 



shocks. Alternatively, high-velocity shocks can result in ion- 
ized lines. 

5.1.1. X-ray emission 

T Tauri stars are known to be strong emitters of moderately 
hard X-rays. Feigelson & Lawson (2004) report X-ray lumi- 
nosities toward young stars in the Chamealeon I north cloud 
of L x = 10 28 - lO^'ergs" 1 whereas Telleschi et al. (2006) find 
X-ray luminosities toward young stars in the Taurus molecu- 
lar clouds of Lx = 10 28 - 10 31 ergs" 1 with a clear stellar mass 
dependence. A number of sources in our sample have been 
identified with X-ray sources and have derived X-ray lumi- 
nosities varying from Lx ~ 4 x 10 28 - 10 31 ergs -1 (Table 4). 

Recently Glassgold et al. (2007) modeled the excitation of 
neon in an X-ray irradiated flaring disk model according to 
D'Alessio et al. (1999) and predict [Nell] and [Nelll] line 
intensities. For a source at 140 pc (the assumed distance 
in Glassgold et al. (2007)) [Nell] line fluxes of (0.6- 1) x 
10" 14 ergcnrV 1 and [Nelll] line fluxes of (5-6.5) x 10" 16 
erg cm -2 s -1 are predicted. The predicted [Ne II] line strength 
of ~ 4-8 x 10~ 6 L Q falls within the observed range of line 
strengths (see Table 4). The predicted [Nelll] line strength 
is at the lower end of the Spitzer IRS [Nelll] upper limits. 
Interestingly the predicted [Nell] line strength agrees well 
at 4000 K (the assumed effective temperature in the Glass- 
gold et al. (2007) model) with the trend in the observed line 
strengths as shown in Figure 9. Also the [Ne II]/[Ne III] line 
ratio for Sz 102 agrees with the line ratios from the Glassgold 
et al. (2007) model. 

The Glassgold et al. (2007) model is based on an assumed 
neon abundance of 10" 4 and an X-ray luminosity and spectral 
temperature of L x = 2 x 10 30 ergs _1 and kT x = 1 keV, appropri- 
ate for solar-mass pre-main sequence stars observed in Orion 
(Wolk et al. 2005). Lower mass, older, and accreting stars 
may have a lower X-ray luminosity (e.g. Feigelson & Law- 
son 2004; Preibisch et al. 2005; Telleschi et al. 2006) leading 
to lower expected line intensities whereas higher X-ray lumi- 
nosities, higher neon abundances (Drake et al. 2005; Cunha 
et al. 2006) or the inclusion of H atom collisions (not included 
in the Glassgold et al. (2007) calculations) may yield higher 
predicted line intensities. The neon line fluxes also scale with 
the disk mass surface density and are therefore sensitive to the 
disk geometry, e.g. flaring or non-flaring disks. 

Of the sources with [Ne II] emission about 30 % are iden- 
tified as X-ray sources. The remaining sources may have es- 
caped detection due to incomplete or sensitivity limited X-ray 
searches or due to source geometry prohibiting the detection 
of the X-rays. A more targeted deep X-ray search would be 
required to confirm a direct relation between observed [Ne II] 
emission and X-ray luminosities. Overall, variations in X-ray 
luminosities, age, stellar mass and geometry appear able to 
cover the two to three orders of magnitude range in the ob- 
served [Ne II] line fluxes. 

5.1.2. EUV radiation 

EUV radiation originating from the stellar chromosphere or 
the accretion shock may be an additional heating component 
of the disk surface and contribute to the neon excitation. EUV 
photons are however quickly absorbed by atomic hydrogen in 
the accretion column (Alexander et al. 2005) and for strong 
accretors will not reach the disk surface. However, for tran- 
sitional objects like T Cha the accretion column can become 
optically thin to EUV photons and some EUV radiation may 



c2d IRS spectra of T Tauri stars: gas-phase lines 



9 



■ — ' 

<N 
I 



cn 
o 



-2 
-4 

-6 
-2 

-4 

-6 
-2 

I -4 



o 



o 



-6 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




i i 1 1 1 1 1 1 1 1 i 1 1 1 1 1 1 i 1 1 1 1 1 1 1 1 i 1 1 1 1 i#r i i i 



1 





0.01 1 
log L IR 



■4 -2 4000 7000 10000 
a T. [K] 



Fig. 9. — Comparison of the observed H2 S(2), [Nell], and [Fel] line strengts in solar luminosities with the mid-IR luminosity, the 12.8-26/im 
spectral index, and the effective temperature of the sources. Plotted with the large black symbols are the line detections, while the small gray 
symbols and arrows are all non-detections. The horizontal dotted lines in the H2 and [Nell] panels show the predicted line strength by Nomura 
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escape the immediate surroundings of the star and reach the 
disk surface, potentially contributing to the ionizing radiation 
at the factor of 2 level. The [Ne II]/[Ne III] ratio is expected to 
increase with the additional EUV contribution given the high 
(41 .0 eV) ionization potential of Ne + . For Sz 102 the tentative 
[Ne III] detection results in a [Ne II]/[Ne III] ratio consistent 
with X-ray excitation. For all other sources the [Ne III] upper 
limits do put any constraints on this (see Fig. 10), however. 

5.1.3. Disk shocks 

The presence of strong [Ne II] could also indicate a possi- 
ble origin of the observed [Ne II] emission induced by shocks 
in the disk. Hartmann & Raymond (1989) describe shocks 
resulting from the stellar wind striking the disk surface at an 
oblique angle. For typical wind velocities of 200 kms" 1 shock 
velocities along the disk surface are estimated to be ~30-40 
kms" 1 . For such shocks and a medium with a density of 



10 5 - 10 6 cm" 3 Hollenbach & McKee (1989) predict [Nell] 
line strengths of <~ 10" 6 - 10~ 4 erg cm" 2 s" 1 sr" 1 . For a 100 AU 
disk at lOOpc, the upper end of this range implies a [Nell] 
line flux of approximately 10~ 14 erg cm" 2 s" 1 , of the same 
order as observed. However, Hollenbach & McKee (1989) 
also predict the H 2 0-0S(1), S(2), and S(3) lines and the 
[Fel]24/mi, [SI] 25.25/im, [FeII]26/xm, and [Si II] 34.8/xm 
spectral lines to be stronger than the [Ne II] line by 1 - 3 orders 
of magnitude for these velocities as illustrated in the lower 
right plot of Fig. 10. The top plots of Fig. 10 show the ob- 
served ratios and lower limits. For almost all of the line pairs, 
the ratios are roughly equal to or larger than unity, except for 
a small number of cases. 

To account for line ratios equal to or larger than unity, 
higher velocity shocks would be required. At high shock ve- 
locities (v > 70kms _1 ) and high densities (> 10 5 cm~ 3 ) the 
J-shock models presented in Hollenbach & McKee (1989) 
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give [Nell], [Fel] 24 /xm, [SI] 25.25 /xm, [Fell] 26 /xm, and 
[Sill] 34.8 /xm lines all of comparable strength of ^0.004 
ergcm" 2 s _1 sr -1 (see bottom right plot in Figure 10). To pro- 
duce a line flux of ~ 10~ 14 ergcm" 2 s _I the shocked emission 
would have to come from a region of the disk with a radius of 
~10 AU at a distance of 100 pc. Higher velocity shocks may 
produce the observed line flux ratios and lower limits, but a 
possible origin for such high velocity shocks is unclear. An- 
other problem lies in the non-detection of the [S I] 25.25/xm, 
the [Fell] 26/xm, and the [Sill] 34.8/xm lines. In particular 
the [S I] line should be detected if high velocity shocks are 
the origin of the observed line emission: the [S I] line is pre- 
dicted to be stronger whereas the detection limits for [Fel] 
and [S I] are comparable in our data (see Fig. 4). 

5.2. [Fel] and [SI] 

As mentioned in Sect. 5.1 the detection of [Fel] in combi- 
nation with the non-detection of other atomic lines, in partic- 
ular [SI] and [Sill], is significant. 

Gorti & Hollenbach (2004) modeled the line emission 
from intermediate aged optically thin disks around G and K 
stars. For disks with low gas masses (10~ 3 to 10~ 2 Mj), the 
[SI] 25.2/xm, [Fell] 26/xm, and [Sill] 35.4/xm lines are ex- 
pected to be the strongest mid-infrared emission lines. How- 



ever as the disk mass increases a larger fraction of the sul- 
fur turns molecular and the [Si II] and [Fe II] emission be- 
comes optically thick. At the same time the [Fe I] lines are 
predicted to become increasingly stronger and will at some 
point, around a O.lMj, dominate over the [S I] emission. Al- 
though specific calculations for these optically thick disks 
are lacking, the seven sources (WXCha, C7-11, HD 132947, 
LkHa271, CokuTau3, Haro 1-17, and XXCha) which show 
strong [Fe I] emission may well have optically thick massive 
gas-rich disks. 

5.3. Molecular hydrogen 

The third significant result of our survey is the non- 
detection of the H 2 0-0 S(0) and S(l) lines for 76 T Tauri and 
Herbig Ae/Be stars. This puts constraints on the mass of warm 
(T - 100-200 K) H 2 gas in the disks around these stars of 
typically a few Mj as illustrated in Figure 8. Models of disk 
heating by stellar UV photons show that the gas temperature 
in the surface layers can be significantly higher than that of 
the dust down to an optical depth for UV photons of ~1. The 
precise temperatures depend on the model details, in particu- 
lar the presence of PAHs, the grain size, the gas/dust ratio, and 
the presence of excess UV over that of the stellar photosphere 
(Jonkheid et al. 2004; Kamp & Dullemond 2004; Nomura & 
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Millar 2005; Nomura et al. 2007; Jonkheid et al. 2006, 2007). 
For interstellar-sized grains (~0.1 /im), the models readily 
give surface temperatures of 100 K or more out to at least 
100 AU. Even models in which the dust grains have grown 
and settled to the midplane have warm surface layers as long 
as some PAHs are still present. The total mass contained in 
this warm layer is however small, ~ 1 % or less of the total 
disk mass. For the specific 0.07 M Q disk studied by Jonkheid 
et al. (2004) the mass at T > 100K is -0.7 M s . Thus, for a 
typical disk mass of 0.01 M©, this may be as low as 0.1 Mj, 
below our upper limits. 

The H2 line fluxes from a protoplanetary disk representa- 
tive of that around TW Hya have been modeled by Nomura 
& Millar (2005). These models include not just thermal exci- 
tation but also UV pumping of the H2 levels. Nomura et al. 
(2007) include X-ray irradiation and the effect of grain size 
distributions. For the higher excitation lines the grain size dis- 
tribution is particularly important with a distribution toward 
smaller grain sizes producing higher line strengths. This re- 
sult is consistent with that of Jonkheid et al. (2007), who find 
lower temperatures in models with grain growth. Figure 7 
shows the distribution of the observed line strengths with the 
predicted line strengths for the Nomura et al. (2007) model in- 
cluded for grain size distributions with maximum grain sizes 



of flmax = 10/um, 1mm, and 10cm. It is seen that the observed 
upper limits are all consistent with this model, even if ex- 
cess UV is included. Note, however, that pure rotational H2 
line fluxes are extremely sensitive to the model details: small 
changes in the heating and cooling processes, as well as the 
treatment of the H/H2 transition zone, can result in signifi- 
cant differences in gas temperatures and an order of magni- 
tude variation in predicted line fluxes (see discussion in Li 
et al. 2002; Roellig et al. 2007). Therefore, comparison of the 
total mass of warm gas between models and observations is 
equally relevant. 

Toward six sources (—8 % of the sample), H2 0-0 S(2) 
and/or S(3) emission is observed, which provides evidence for 
the presence of a significant hot (T > 500 K) gas component 
in the disks. Hot gas (T > 500 K) is observed toward a num- 
ber of sources through the H 2 0-0 S(2) line, most convincingly 
toward Sz 102, EC 82, CedIRSIRS6, EC 74, and EC 92, For 
Sz 102 a number of higher transition lines are also seen. The 
observed H2 0-0 S(2) line strengths are more than a factor of 
10 higher than those predicted in Nomura & Millar (2005); 
Nomura et al. (2007) (see Fig. 7). Given the non-detection of 
the S(l) line toward the same sources, which is predicted to 
have a similar strength, this is an indication that these disks 
have an additional source of emission from hot molecular hy- 
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drogen. The upper limits for H2 0-0 S(2) and S(3) are higher 
for almost all sources than the Nomura & Millar (2005) pre- 
dicted values. An additional hot component may therefore be 
present in these sources as well below our detection limit. 

None of the sources with evidence for an additional hot 
component show evidence for PAH emission (see Geers et al. 
2006). Of the eight sources with detected H 2 0-0S(2), four 
show strong [Nell] emission, giving support to the idea of 
a common heating and excitation mechanism through X-rays 
or EUV. Of the [Ne II] sources two, T Cha and RR Tau, show 
strong PAH emission. Considering the limitations in observ- 
ing the PAH emission as described by Geers et al. (2006) we 
can at this stage draw no conclusions about the relation be- 
tween the hot H2 emission, [Ne II] emission, and the impor- 
tance of PAHs and small grains. 

An origin in a high velocity shock as discussed in Section 
5.1.3 could produce the enhanced S(2) and S(3) line strengths 
while keeping the line strengths of S(0) and S(l) reduced. 
However, as discussed in Section 5.1.3, the main problem 
with invoking shocks to explain [Ne II] and H2 is to accom- 
modate both the detections of [Fe I] and the non-detections 
of [S I] and [Fe II] 26/im. Figure 1 1 shows the observed and 
model line ratios with respect to H2, similar those in Fig. 10. 
An origin in an oblique stellar wind shock faces the same 
problems as discussed for [Nell] in Section 5.1.3. 

6. CONCLUSIONS 

A survey of the mid-infrared gas phase pure rotational lines 
of molecular hydrogen and a number of atomic fine structure 
transitions has been carried out toward a significant sample 
of 76 circumstellar disks with the Spitzer IRS. The principal 
findings include: 

• [Ne II] is detected toward ~ 20 % of the sources and 
[Ne III] tentatively in one source. The [Ne II] detections 
and the [Ne II]/[Ne III] line flux ratio are consistent with 
disk heating and excitation of [Ne II] through X-rays as 
presented in Glassgold et al. (2007). Excitation through 
EUV radiation may contribute. Better constraints on 
the X-ray luminosities and [Ne III] fluxes are required 
to distinguish the two contributions. 

• [Fe I] is detected toward ~ 9 % of the sources. No other 
low excitation atomic lines, such as [Fe II], [SI] and 
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[Sill], are detected. This suggests that these sources 
may possess optically thick disks with gas masses of at 
leastaO.lMj. 

• Except for a tentative detection toward Sz 102, no com- 
pact H2 0-0S(0) and S(l) emission is observed toward 
any of the sources in our sample, setting limits of a few 
Jovian masses on the mass of the warm r ex = 100 K gas 
in the disks. These limits are above model predictions. 
The H 2 line flux upper limits are also consistent with 
recent T Tauri disk model predictions by for example 
Nomura et al. (2007). Earlier tentative ISO detections 
of H2 in two Herbig Ae disks are not confirmed. 

• Hot (T > 500 K) H2 gas has been detected toward ~ 
8 % of the sources. Given the high upper limits for the 
rest of the sources, the fraction may be higher. The 
detection of the hot gas suggests the presence of an ad- 
ditional source of hot H2 emission not included in the 
most recent disk models (e.g. Nomura et al. 2007). 

• An origin of the enhanced H2 emission in oblique 
shocks due to winds interacting with the disk surface is 
not consistent with the non-detection of atomic lines, in 
particular the non-detections of [S I] and [Fe II] 26^m. 

The bright [Nell] lines detected at 12.8 /im are excel- 
lent targets for follow-up observations with high dispersion 
echelle spectrometers on 8-10m class telescopes (TEXES, 
VISIR). The measured spatial profiles and line shapes would 
provide exacting tests of the X-ray mediated disk emission 
proposed here and could definitely rule out any high-velocity 
shock mechanism. 
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178 


1.9-40.0 


M1-M2 


2,12,13,14,15,17 


30 


HD 98922 


ll h 22 m 31 s .7 


-53° 22' 11' 


'.3 


5640704 


(-/6*1/14*1)*2 


HAeBe 


>540 


- 


B9 


2 


31 


HD 101412 


ll h 39 m 44 s .5 


-60° 10' 27' 


'.8 


5640960 


(14*1/31*1/60*1)*2 


HAeBe 


160 


- 


B9.5 


2 


32 


TCha 


ll h 57 ra 13 s .5 


-79° 21' 31' 


'.3 


5641216 


(-/31*2/60*1)*2 


TTs 


66 


4.1-12.5 


G2-K0 


2,3,12,13,15 


33 


IRAS 12535-7623 


12 h 57 m ll s .8 


-76° 40' 11' 


'.6 


11827456 


(-/121*l/60*2)*2 


TTs 


178 


0.2 


MO 


2,13,15 


34 


Sz50 


13 h 00 m 55 s .4 


-77° 10' 22' 


'.1 


11827456 


(-/121*l/60*2)*2 


TTs 


178 


1.7 


K7-M3 


2,13,15 


35 


ISO-Chall 54 


13 h 00 m 59 s .2 


-77° 14' 02' 


'.8 


15735040 


( 14*2/3 1*4/-)*2 


BD 


178 


- 




21 


36 


DL Cha 


13 h 06 m 08 s .4 


-77° 06' 27' 


'.4 


5642240 


(-/6*2/6*5)*2 


TTs 


178 


- 


M6 


2,13,15,22 


37 


HD 132947 


15 h 04 m 56 s .O 


-63° 07' 52' 


'.7 


5643008 


(14*1/31*2/241*2)*2 


HAeBe 


>60 


- 


AO 


2,13,15,22 


38 


HD 135344 


15 h 15 ra 48 s .4 


-37° 09' 15' 


'.8 


5657088 


(-/121*l/60*2)*2 


HAeBe 


140 


17.0 


A0-F4 


2 


39 


HT Lup 


15 h 45 m 12 s .9 


-34° 17' 30' 


'.5 


5643264 


(14*1/31*1/60*1)*2 


TTs 


145 


0.4-2.3 


K2 


2,3,14,23,24 


40 


HT Lup 


15 h 45 m 12 s . 9 


-34° 17' 30' 


'.5 


9829120 


(14*2/31*2/60*2)*1 


TTs 


145 


0.4-2.3 


K2 


2,3,14,23,24 


41 


GWLup 


15 h 46 m 44 s .7 


-34° 30' 35' 


'.3 


5643520 


(14*l/121*2/60*4)*2 


TTs 


100 


1.3-3.4 


M2-M4 


2,3,24 


42 


Sz73 


15 h 47 m 57 s .O 


-35° 14' 35' 


'.2 


5644032 


(14*1/31*1/60*1)*2 


TTs 


100 


2.6-5.4 


K2-Ma 


2,3,24 


43 


GQ Lup 


15 h 49 m 12M 


-35° 39' 05' 


'.() 


5644032 


(14*1/31*1/60*1)*2 


TTs 


100 


0.1-1.0 


K7-M0 


2,3,14,23,24 


44 


IMLup 


15 h 56 m 09 s .2 


-37° 56' 06' 


'.4 


5644800 


(14*1/31*2/60*1)*2 


TTs 


140 


0.1-1.5 


MO 


2,3,20,24 


45 


RULup 


15 h 56 m 42 s .3 


-37° 49' 15' 


'.6 


5644800 


(14*1/31*2/60*1)*2 


CTTs 


140 


0.1-0.8 


K3-M0 


2,23,24 


46 


RYLup 


15 h 59 m 28 s .4 


-40° 21' 51' 


'.1 


5644544 


(14*1/31*1/60*1)*2 


TTs 


150 


1.6-3.2 


K0-K4 


2,3,14,24 


47 


EX Lup 


16 h 03 ra 05 s .5 


-40° 18' 24' 


'.8 


5645056 


(14*1/31*1/60*1)*2 


TTs 


150 


1.4-5.6 


MO 


2,3,14,24 



TAB LEI — Continued 



if 


Source 






rtUK 


• 

tint [sj 




u 


Age 


Spectral type 


b 

Kererences 












(SL/SH/LH)*ndith. 




TdcI 


TMvrl 






[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


48 


Szl02 


16 h 08 m 29 s .7 


-39° 03' 11 ".2 


9407488 


(14*1/121*2/60*1)*2 


TTs 


200 


- 


M0 


2,3,24 


49 


AS 205 


16 h ll ra 3P.3 


-18°38'26".2 


5646080 


(6*2/6*2/14*2)*2 


TTs 


120 


0.1-0.4 


K5 


2,3,20,25 


50 


Haro 1-1 


16 h 21 m 34 s .7 


-26° 12' 27".0 


9833472 


(14*1/31*4/60*1)*2 


CTTs 


125 


13.1 


K5-K7 


3,26,27,28 


51 


Haro 1-4 


16 h 25 m 10 s .5 


-23° 19' 14". 5 


9833216 


(-/31*2/60*1)*2 


TTs 


125 


0.8 


K4-K6 


26 


52 


DoAr 24E 


16 h 26 m 23 s .4 


-24° 21' 00" .0 


5647616 


(-/31*1/14*2)*2 


TTs 


125 


1.5-2.5 


K0-K1 


2,20,25,29,30,31,32 


53 


SR21 


16 h 27 m 10 s .3 


-24° 19' 12" .4 


5647616 


(-/31*1/14*2)*2 


TTs 


125 


1.0-3.0 


F4-G2.5 


2,25,30,31,33,34,35 


54 


IRS51 


16 h 27 m 39 s .8 


-24° 43' 15". 2 


9829888 


(-/31*2/14*2)*2 


2 


125 


0.8 


G5-K7 


3,16,25,33 


55 


SR9 


16 h 27 m 40 s .3 


-24° 22'04".l 


12027392 


(-/31*l/60*2)*2 


TTs 


125 


0.8 


K5-M2 


2,3,25,26,29,33,34 


56 


V8530ph 


16 h 28 m 45 s .3 


-24° 28' 18". 8 


12408576 


(-/31*4/60*4)*2 


TTs 


125 


0.8 


M1.5 


25,26 


57 


ROX42C 


16 h 31 m 15 s .7 


-24° 34' 01". 9 


6369792 


(-/31*2/60*1)*2 


TTs 


125 


3.7 


K4-K6 


16,25 


58 


ROX43A 


16 h 31 m 20 s .l 


-24° 30'05".0 


15914496 


(-/31*l/60*2)*2 


TTs 


125 


12.2 


GO 


25,34,36 


59 


IRS60 


16 h 31 m 30 s .9 


-24° 24' 39". 6 


6370048 


(-/31*1/60*1)*2 


- 


125 


- 




16,25,36,37 


60 


Haro 1-16 


16 h 31 m 33 s .5 


-24° 27'37".l 


12664064 


(-/31*1/60*1)*2 


TTs 


125 


0.5-1.9 


K2-K3 


14,20,25,36 


61 


Haro 1-17 


16 h 32 m 21 s .9 


-24° 42' 14". 8 


11827712 


(60*l/121*4/60*4)*2 


TTs 


125 


- 


M2.5 


2,3,25 


62 


RNO 90 


16 h 34 m 09 s .2 


-15° 48' 16". 9 


5650432 


(14*1/31*1/14*2)*2 


TTs 


140 


6.0 


G5 


1,2,23,25,38,39 


63 


WaOph6 


16 h 48 m 45 s .6 


-14° 16' 35". 8 


5650688 


(14*1/31*1/60*1)*2 


2 


140 


0.9 


K7 


25,37 


64 


VI 121 Oph 


16 h 49 ra 15 s .3 


-14° 22' 08".8 


5650688 


(14*1/31*1/60*1)*2 


CTTs 


125 


0.5 


K5 


3,23,25,37 


65 


HD 163296 


17 h 56 m 21 s .3 


-21° 57' 22" .0 


5650944 


(-/6*2/14*2)*2 


HAeBe 


122 


4.0-6.0 


A0-A2 


2 


66 


HD 163296 


17 h 56 m 21 s .3 


-21° 57' 22" .0 


9830144 


(-/6*4/14*4)*l 


HAeBe 


122 


4.0-6.0 


A0-A2 


2 


67 


VV Ser 


18 h 28 m 47 s .9 


0° 08' 39". 8 


5651200 


(6*2/31*l/60*l)*2 


HAeBe 


260 


0.6 


B1-A3 


2 


68 


SSTc2dJ182900.9+02931 


18 h 29 m 00 s .9 


0° 29'31".6 


13210112 


(14*l/31*2/-)*2 


- 


260 


- 




- 


69 


SSTc2d Jl 82909.8+03446 


18 h 29 m 09 s .8 


0° 34' 45". 8 


13210624 


(14*l/31*l/-)*2 


- 


260 


- 




- 


70 


SSTc2dJl 82928.2+02257 


18 h 29 m 28 s .2 


-0° 22' 57" .4 


13210368 


(14*l/31*2/-)*2 




260 








71 


EC74 


18 h 29 m 55 s .7 


1° 14'31".6 


9407232 


(14*l/121*l/60*2)*2 


2 


260 






26 


72 


EC82 


18 h 29 m 56 s .9 


1° 14' 46". 7 


9407232 


(14*l/121*l/60*2)*2 


TTs 


260 




M0 


2,4,26,40 


73 


EC90 


18 h 29 ra 57 s .7 


1° 14'06".0 


9828352 


(6*2/31*l/14*l)*2 


TTs 


260 






4,26,40,41 


74 


EC92 


18 h 29 m 57 s .9 


1° 12'51".5 


9407232 


(14*l/121*l/60*2)*2 


TTs 


260 




K7-M2 


18,41,42,43 


75 


CK4 


18 h 29 m 58 s .2 


1° 15'21".6 


9407232 


(14*l/121*l/60*2)*2 


TTs 


260 


6.8 


K3 


2,26,40,43,44 


76 


LkHct 348 


18 h 34 ra 12 s .6 


-0° 26' 21". 8 


9831424 


(6*2/6*2/14*2)*2 




260 









a Assumed cloud distances; Chamaeleon (178 pc Whittet et al. 1997) Lupus I, IV (150 pc) and Lupus III (200 pc) Comeron in prep., Ophiuchus (125 pc de Geus et al. 1989), Perseus (250 pc, see discussion in 
Enoch, Young, Glenn, Evans, Golwala, Sargent, Harvey, Aguirre, Goldin, Haig, Huard, Lange, Laurent, Maloney, Mauskopf, Rossinot, & Sayers 2006), Taurus-Auriga (160 pc Kenyon et al. 1994), Serpens (260 pc 
Straizys, Cernis, & Bartasiute 1996), except for T Cha, HD 98922, & HD 163296 (van den Ancker et al. 1998), BF Ori, (Parenago 1954), HD 101412 & HD 135344 (de Zeeuw et al. 1999), & HD 132947 (minimum 
Tycho distance) 

b A non-exhaustive literature reference list: (1) Clark (1991), (2) Kessler-Silacci et al. (2006), (3) Herbig & Bell (1988), (4) Casali & Eiroa (1996), (5) Ladd et al. (1993), (6) Beichman et al. (1984), (7) Myers 
et al. (1987), (8) Kenyon et al. (1998), (9) Edwards et al. (2006), (10) Beckwith et al. (1990), (11) Sargent et al. (2006), (12) Lawson et al. (1996), (13) Gauvin & Strom (1992), (14) Appenzeller et al. (1983), (15) 
Hartigan (1993), (16) Beichman et al. (1988), (17) Cambresy et al. (1998), (18) Pollanen & Feldman (1995), (19) Persi et al. (2000), (20) Schegerer et al. (2006), (21) Allers et al. (2006), (22) Lopez & Girard 
(1990), (23) Carballo et al. (1992), (24) Hughes et al. (1994), (25) Andre & Montmerle (1994), (26) Churchwell & Koornneef (1986), (27) Rydgren (1980), (28) Meyer et al. (1993), (29) Elias (1978), (30) Lada & 
Wilking (1984), (31) Ward-Thompson (1993), (32) Barsony et al. (1997), (33) Wilking et al. (1989), (34) Jensen & Mathieu (1997), (35) Brown et al. (in prep.), (36) Greene et al. (1994), (37) Ichikawa & Nishida 
(1989), (38) Reipurth et al. (1996), (39) Parker (1991), (40) Kaas (1999), (41) Haisch et al. (2002), (42) McMullin et al. (2000), (43) Hurt & Barsony (1996), (44) Zhang et al. (1988). 



TABLE 2 

Observed linefluxes and 1<t uncertainties (l(T' 6 ergcrrr 2 s~' ) 



# 


Source 


Dis. [pc] 


H 2 S(0) 


H 2 S(l) 


H2 S(2) 


H 2 S(3) 


[Ne II] 


[Ne III] 


[Fe I] 


[Fe II] 1 8 


[Fe II] 25 


[S I] 


[Si II] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


1 


RNO 15 


250 


- (48) 


- (63) 


-(30) 


- (470) 


- (85) 


- (69) 

yys / 


- (84) 


- (78) 


- (44) 

WW 


- (43) 


-(43) 


2 


LkHo 270 


250 


- (35) 


- (42) 


-(27) 


-(63) 


-(29) 


-(21) 


-(70) 


- (27) 


- (28) 


- (35) 


- (180) 


3 


LkHa 271 


250 


- (53) 


- (10) 


-(12) 


-(42) 


-(12) 


-(8) 


880 (130) 


- (16) 


- (19) 


- (23) 


-(150) 


4 


LkHa 326 


250 


- (39) 


- (25) 

WW 


- (22) 

WW 


- (110) 


- (21) 


- (26) 


- (88) 

WW 


- (24) 


- (15) 

WW 


- (33) 


- (160) 


5 


LkHo 327 


250 


-(63) 


-(60) 


-(53) 


- (240) 


-(64) 


- (48) 


-(76) 


-(59) 


-(34) 


- (52) 


-(48) 


6 


LkHo 330 


250 


- (110) 


- (45) 


- (27) 

W ' / 


- (220) 


- (18) 

W W 


- (34) 


- (160) 


- (42) 
ww 


- (100) 

W "W 


- (100) 
V 1 w w 


- (89) 

WW 


7 


IRAS 03446+3254 


250 


-(76) 


-(10) 


-(6) 


-(62) 


36 (6) 


-(13) 


-(110) 


-(13) 


-(44) 


-(46) 


-(110) 


8 


LkCa 8 


140 


-(22) 


-(15) 


-(11) 




-(10) 


-(9) 


-(71) 


-(16) 


-(19) 


-(27) 


-(97) 


9 


IQ Tau 


140 


- (38) 

WW 


- (27) 

W ' / 


- (24) 

WW 




- (25) 


- (23) 

WW 


-(61) 


- (20) 

WW 


- (18) 
w w 


- (28) 

WW 


- (140) 
V *^w 


10 


FX Tau 


140 


- (37) 


-(24) 


-(34) 




-(19) 


-(17) 


-(97) 


-(29) 


-(25) 


- (19) 


- (130) 


1 1 


V710Tau 


140 


- (66) 


- (12) 


- (12) 




-(8) 


-(11) 


- (150) 


- (11) 
y 1 1 j 


- (28) 


- (25) 


- (72) 


12 


DN Tau 


140 


- (33) 


- (18) 


- (21) 




-(14) 


- (17) 


- (77) 


- (18) 


- (23) 


- (31) 


- (83) 


13 


CoKu Tau 3 


140 


- (29) 

y^s / 


- (19) 


- (14) 




- (21) 


-(16) 


450 (80) 


- (23) 


- (26) 

WW 


- (18) 


- (62) 

WW 


14 


CoKu Tau 4 


140 


- (49) 


- (30) 


- (24) 




- (10) 


-(13) 


-(75) 


- (33) 

WW 


- (33) 

WW 


- (36) 

WW 


- (87) 

W ' / 


15 


BF Ori 


400 


- (92) 


- (38) 


- (48) 


- (1000) 


- (48) 

WW 


- (44) 

WW 


- (73) 


- (49) 


- (33) 

WW 


- (180) 

W U W 


- (51) 


16 


RR Tau 


160 


- (81) 


- (52) 


- (48) 


- (790) 


- (63) 


- (46) 


-(77) 


- (52) 

WW 


- (67) 

W ' 1 


- (67) 

W ' / 


-(73) 


17 


IRAS 08267-3336 


400? 


- (35) 


- (20) 


- (22) 


- (82) 


77 (22) 
' ' ww 


-(22) 


- (73) 


- (32) 

WW 


- (26) 

WW 


- (23) 

WW 


- (32) 


18 


SX Cha 


178 


- (46) 


- (28) 


- (19) 




- (27) 

W ' > 


- (20) 


- (54) 


- (27) 

W ' / 


- (24) 

WW 


- (30) 

WW 


- (54) 


19 


SY Cha 


178 


- (36) 


- (13) 


- (12) 




- (12) 
w w 


-(11) 


- (59) 

w w 


- (21) 
w w 


- (35) 

WW 


- (19) 

w w 


- (44) 


20 


TW Cha 


178 


- (53) 


- (18) 


- (14) 




- (ii) 


- (16) 


-(79) 


- (15) 


- (27) 


- (34) 


-(50) 


21 


Ced 110 IRS6 


178 


- (66) 

WW 


- (26) 


49 (11) 




50 (13) 


-(21) 


- (60) 

WW 


- (33) 

WW 


- (52) 

WW 


- (47) 


-(48) 

WW 


22 


B35 


178 


-(49) 


- (11) 


-(8) 




-(4) 


-(14) 


-(56) 


- (10) 


- (17) 


-(23) 


-(48) 


23 


VW Cha 


178 


-(65) 


- (66) 


-(57) 




280 (38) 


-(55) 


-(99) 


-(49) 


-(30) 


- (87) 


-(80) 


24 


VZ Cha 


178 


- (53) 

WW 


- (23) 


- (35) 




- (26) 

WW 


- (30) 

WW 


- (72) 


- (22) 


- (30) 

WW 


- (27) 


- (43) 

WW 


25 


WX Cha 


178 


-(27) 


-(28) 


-(22) 




-(16) 


-(26) 


580 (66) 


-(23) 


-(46) 


-(24) 


-(41) 


26 


ISO-Cha237 


178 


- (36) 


- (19) 


33 (9) 

• J -' w/ 




- (20) 

WW 


- (15) 

WW 


-(87) 

W ' t 


- (27) 

W ' / 


- (31) 
w w 


- (33) 

WW 


- (46) 

WW 


27 


C7-11 


178 


-(21) 


-(11) 


-(8) 




-(6) 


-(11) 


470 (54) 


-(13) 


-(32) 


-(16) 


-(37) 


28 


HM 27 


178 


- (34) 

WW 


- (17) 


- (12) 




- (10) 


- (21) 


- (64) 

WW 


- (17) 
y 1 1 / 


- (43) 

WW 


- (21) 
w w 


- (40) 

WW 


29 


XX Cha 


178 


- (34) 


- (9) 


- (7) 




35 (7) 


-(9) 


340 (81) 


- (12) 


- (22) 


- (20) 

WW 


- (52) 

WW 


30 


HD 98922 


> 540 


- (350) 


- (850) 


- (840) 




- (1400) 


- (750) 


- (600) 


- (830) 


- (370) 


- (450) 


- (190) 


31 


HD 101412 


160 


- (74) 


- (89) 


- (90) 

WW 


- (690) 

W-' w 


- (120) 


- (74) 


- (100) 


- (110) 


- (32) 


- (65) 


-(66) 


32 


T Cha 


66 


- (47) 


- (17) 


-(22) 




320 (21) 


- (20) 


- (60) 


- (29) 


- (28) 


- (44) 


-(50) 


33 


IRAS 12535-7623 


178 


- (53) 


- (15) 


- (12) 




-(13) 


- (13) 


-(67) 


- (18) 


- (23) 

WW 


- (25) 

WW 


- (46) 

WW 


34 


Sz5() 


178 


- (43) 


- (17) 


- (14) 




- (13) 


-(13) 


- (45) 

WW 


- (14) 


- (40) 

WW 


- (21) 
w w 


- (70) 

V ' w 


35 


ISO-Chall 54 


178 




- (12) 


- (16) 


W ' > 


- (13) 


- (11) 




- (24) 








36 


DL Cha 


178 


- (300) 


- (1000) 


- (1600) 




- (2100) 


- (1200) 


- (410) 


- (1100) 


- (400) 

W W W 


- (380) 


- (140) 
w^w 


37 


HD 132947 


> 60 


- (2) 


- (7) 


-(17) 


-(86) 


- (11) 


-(19) 


45 (6) 


- (21) 


-(5) 


-(3) 


-(18) 


38 


HD 135344 


140 


- (150) 


-(51) 


-(46) 




-(32) 


-(35) 


- (160) 


- (67) 


- (110) 


- (120) 


-(160) 


39 


HT Lup 


145 


- (80) 

WW 


- (110) 


- (47) 


- (600) 


- (96) 
w w 


- (75) 


- (130) 


- (81) 


- (78) 


- (91) 


- (74) 

V ' w 


40 


HT Lup 


145 


- (160) 


- (220) 


- (130) 


- (1100) 


- (120) 


- (180) 


- (210) 


- (240) 


- (140) 


- (130) 


-(110) 


41 


GW Lup 


100 


- (42) 

WW 


- (10) 


-(6) 

VW 


- (62) 

WW 


- (7) 


- (14) 


- (39) 

WW 


- (27) 
y^ 1 j 


- (48) 

WW 


- (13) 


-(31) 


42 


Sz73 


100 


-(47) 


-(33) 


-(35) 


-(150) 


160 (24) 


-(35) 


-(63) 


-(36) 


-(37) 


-(41) 


-(49) 


43 


GQ Lup 


100 


-(31) 


-(36) 


-(38) 


-(180) 


-(27) 


-(38) 


-(67) 


-(46) 


-(42) 


-(30) 


-(69) 


44 


IM Lup 


140 


- (45) 


- (32) 

WW 


- (20) 

WW 


- (220) 


61 (16) 


-(27) 

W ' / 


- (55) 

WW 


- (31) 

W W 


- (34) 

WW 


- (23) 

WW 


-(36) 

WW 


45 


RULup 


140 


-(110) 


- (100) 


- (130) 


- (720) 


-(92) 


- (120) 


-(120) 


- (120) 


-(93) 


-(150) 


-(73) 


46 


RYLup 


150 


-(88) 


-(55) 


-(40) 


- (370) 


-(32) 


-(42) 


-(110) 


-(HO) 


-(82) 


-(73) 


-(89) 


47 


EX Lup 


150 


-(52) 


-(38) 


-(33) 


- (230) 


-(22) 


-(40) 


-(88) 


-(37) 


-(45) 


-(36) 


-(100) 


48 


Szl02 


200 


-(35) 


32 (10) 


92 (6) 


690 (56) 


360 (12) 


23 (7) 


-(49) 


49 (13) 


-(9) 


-(9) 


-(25) 


49 


AS 205 


120 


- (250) 


- (350) 


- (340) 


- (2600) 


- (270) 


- (400) 


- (430) 


- (390) 


- (290) 


- (360) 


-(150) 


50 


Haro 1-1 


125 


-(26) 


-(24) 


-(18) 


-(55) 


-(17) 


-(14) 


-(46) 


-(22) 


-(26) 


-(39) 


-(52) 


51 


Haro 1-4 


125 


-(28) 


-(35) 


-(35) 




-(21) 


-(30) 


-(53) 


-(37) 


-(34) 


-(35) 


-(61) 


52 


DoAr 24E 


125 


-(91) 


- (140) 


-(130) 




- (120) 


-(150) 


- (160) 


-(180) 


-(HO) 


-(92) 


- (100) 


53 


SR21 


125 


- (420) 


- (220) 


-(61) 




- (100) 


-(190) 


- (530) 


- (390) 


- (410) 


- (410) 


-(270) 


54 


IRS51 


125 


-(110) 


-(65) 


-(30) 




-(47) 


-(36) 


-(92) 


-(61) 


-(96) 


-(95) 


-(91) 



TABLE 2 — Continued 



# 


Source 


Dis. [pc] 


H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 


[Ne II] 


[Ne III] 


[Fe I] 


[Fe II] 18 


[Fe II] 25 


[S I] 


[Si II] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


55 


SR 9 


125 


- (60) 


-(57) 


-(42) 




-(21) 


-(32) 


-(68) 


-(54) 


-(38) 


-(35) 


-(53) 


56 


V853upn 


125 


-(27) 


-(26) 


-(19) 




1 10 (21) 


-(24) 


-(54) 


-(30) 


-(36) 


-(43) 


-(63) 


57 


ROX42C 


125 


-(37) 


-(31) 


-(23) 




-(18) 


-(38) 


-(41) 


-(28) 


-(20) 


-(39) 


-(79) 


58 


ROX43A 


125 


-(72) 


-(110) 


-(74) 




-(73) 


-(91) 


-(110) 


-(130) 


-(47) 


-(86) 


-(70) 


59 


IRS 60 


125 


-(30) 


-(47) 


-(86) 




141) (4z) 


-(65) 


-(62) 


-(50) 


-(24) 


-(75) 


-(90) 


60 


Haro 1-16 


1 1C 

123 


i A Q\ 

- (48) 


- (67) 


- (36) 




-(35) 


- (49) 


-(110) 


- (57) 


- (24) 


- (73) 


- (83) 


61 


Haro 1-17 


125 


-(33) 


-(13) 


-(6) 


-(59) 


16 (4) 


-(6) 


400 (88) 


-(14) 


-(29) 


-(20) 


-(51) 


62 


RNO 90 


140 


- (140) 


- (190) 


-(130) 


- (730) 


-(94) 


-(160) 


- (170) 


-(160) 


-(120) 


-(160) 


-(110) 


63 


Wa Oph 6 


140 


-(30) 


-(72) 


-(56) 


- (300) 


-(43) 


-(54) 


-(37) 


-(47) 


-(16) 


-(23) 


-(21) 


64 


V 1 1 2 1 (Jpn 


125 


- (56) 


-(170) 


- (120) 


- (1600) 


- (140) 


- (140) 


-(89) 


- (250) 


- (54) 


- (81) 


-(51) 


65 


HD loiivo 


122 


- (340) 


- (610) 


— (460) 




- (440) 


- (370) 


- (450) 


- (720) 


- (260) 


— (400) 


- (170) 


66 


HD 163296 


122 


- (510) 


— (1200) 


— (1 100) 




- (400) 


- (1200) 


- (1200) 


- (1 100) 


- (270) 


- (370) 


- (350) 


67 


VV Ser 


260 


- (100) 


-(190) 


-(170) 


- (2000) 


-(130) 


-(160) 


- (120) 


-(130) 


-(71) 


-(76) 


-(83) 


68 


SS lc2d J 1 82900.9 


260 




-(32) 


-(28) 


-(130) 


-(16) 


-(31) 




-(37) 








oy 


CCT^1^ T1 OTGnO Q 

icza j lozyuy.o 


ZOU 




{A Q\ 

- (4») 


- (43) 


- (230) 


- (44) 


i"i^\ 

-(33) 




- (63) 








70 


SSTc2dJ182928.2 


260 




-(160) 


-(110) 


- (1400) 


800 (120) 


-(150) 




- (260) 








71 


EC74 


260 


-(73) 


-(12) 


46 (7) 


440 (57) 


25 (6) 


-(8) 


-(63) 


-(12) 


-(160) 


-(150) 


-(130) 


72 


EC82 


260 


-(93) 


- (100) 


170 (35) 


-(740) 


-(63) 


-(51) 


-(150) 


-(130) 


-(110) 


-(85) 


-(120) 


73 


EC90 


260 


- (210) 


- (410) 


- (170) 


- (5800) 


- (200) 


-(310) 


-(310) 


- (610) 


-(250) 


- (280) 


-(190) 


74 


EC92 


260 


-(90) 


-(56) 


69 (18) 


- (420) 


100 (18) 


-(38) 


-(89) 


-(54) 


-(48) 


-(76) 


- (160) 


75 


CK4 


260 


-(38) 


-(27) 


-(19) 


-(81) 


-(12) 


-(23) 


-(50) 


-(24) 


-(18) 


-(21) 


-(58) 


76 


LkHa 348 


260 


-(48) 


-(190) 


- (300) 


- (2200) 


- (220) 


-(210) 


-(84) 


- (560) 


-(38) 


-(41) 


-(30) 



TABLE 3 

Diagnostic parameters of the disk gas contents. 



# a 


Source 




■'Vwarm 

(r = 100 AU) L 


A/T b 
'''warm 


7hot b 


A^hot 

/ t*\ AT T\C 

(r = 2 AU) 


^hot 






[K] 


rl/\99 — 9n 

[10 Z2 cm z ] 


[Mj] 


[K] 


„/i?9 — 9n 

[10 2Z cm z ] 


[Kr J Mj] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


I 


RNO 15 




< 74 


< 9.2 




< 75 


< 3.7 


2 


T k Ffrv 770 




< 59 


< 7 3 




< 18 


< 0.9 


1 
j 


T V Urv 971 
i_jK net z / 1 




/ ?.< 
<^ jj 


■s 4 1 




\ lO 


< 8 


4 


T V Ffrv 196 

IjK net jzo 




< 44 


< 5 5 




<^ 17 


<f 1 Q 
<^ l .y 


< 

j 


T V Ffz-v 197 
i_jK net jz / 




< 86 


< 10 6 




< 80 










< 1 10 


< 13.0 




< 72 


< 3.5 


7 


TRAS 03446+^254 




< 43 


< 5.3 




< 30 


< 1.5 


g 






<r 8 






<^ 4 


< 0.2 


9 


IQTau 




<r 1 5 


<r 1 8 






< 0.5 


in 

1 U 


rJ\. Idll 






<r- 1 6 




<v ij 


<r- 7 


1 1 


V7 1 Tan 
V / 1U IdU 




<^ 14 


1 7 




<^ 4 


< 0.2 


12 


DNTau 




<^ 1 1 

X II 


<r 1 4 




<^ 8 




1 3 


rVVKn Tan 1 




/ 1 1 


^ i 9 
<v i . j 




/ < 

<^ J 


/HI 

U.J 


14 


fnKii Tan 4 




<^ 18 


< 2.2 




\ y 


< 0.5 


15 


BF Ori 




< 220 


< 26 8 




< 225 


< 11.2 


I 6 

1 D 


RR Tau 




<^ JZ 






< 155 


/ 7 7 


1 7 


TR A <5 0,8967 1116 




,-• no 


<v 1Z.Z 




nn 

<s / / 


/in 
<_ j.y 


1 8 


oA v^lld 




/ 97 
<^ z / 


1 4 




<^ 1 9 


/ n 6 


19 


SYCha 




< 16 


< 2.0 




< 7 


< 0.4 


20 


TWCha 




< 24 


< 2.9 




< 9 


< 0.4 


21 


fVrl 1 1 TRS6 
LcU 1 1U IJxijO 




<r" 19 
<^ jZ 


/TO 

<^ 




^ in 

S 1U 


< 5 


zz 


R1S 




1 8 

"v, 15 


<. Z.Z 




^ «: 
j 


^ n i 

<. U.J 


91 

Z J 


VW Phi 
V VV V^lld 




/ ,1Q 
<^ to 


<^ J .y 




^ -27 


/ 1 Q 


24 


VZ Cha 




< 26 


<r 1 9 

\ J.Z 




< 22 


<r 1 1 
\ i.i 


Z J 


vv a L.iia 




<r- 90 
zu 


<^ Z.J 




^ i c 

<~. ij 


/ n 7 

<, u. / 


26 


TSO-rha917 




< 20 


<r 9 5 

\ Z.J 




7 


\ u.t 


27 


C7-11 




<r 19 

\ 1Z 


<r 1 4 




\ j 


< 0.3 


28 


HM 27 




< 18 


< 2.3 




< 7 


< 0.4 


29 


XX Cha 




< 14 


< 1.7 




< 4 


< 0.2 


30 


T4D Q8Q99 




< 2200 


< 274.3 




< 11000 


< 546.2 


31 


HD 101412 




< 40 


< 4.9 




< 110 


< 5.4 


32 


TCha 




<r 9 
\ z 


<r 4 




<- 1 
\ i 


<r n 1 


33 


IRAS 12535-7623 




< 22 


< 2.7 




< 7 


< 0.4 


34 


Sz50 




< 21 


< 2.5 




< 8 


< 0.4 


35 


TSO-ChaTI 54 




< 6 


< 0.8 




< 14 


< 0.7 


36 


DT Cha 

17L- V^IlLl 




< 200 


< 24.3 




< 1975 


< 98.5 


37 


T4D 1 19Q47 




< 0.2 


< 0.0 






< 1 


Ik 

JO 


un 1 1S144 




/ A\ 
<^ 4- 1 


<r" S 1 
<^ J . 1 




\ lO 


< 9 


39 


HT Lup 




< 42 


< 5.1 




< 47 


< 2.4 


40 


HT Lup 




< 84 


< 10.4 




< 75 


< 3.8 


41 


nWT un 

V 1 VV LUU 




<" 4 


< 0.6 




< 4 


< 0.2 


42 






/ Q 


<r" 1 1 




v. o 


/ 04 


41 








< . VJ.7 




/ Q 
<^ o 


/ 04 


44 


TMT nn 




<r 17 


<r 9 1 

<^ Z. 1 




<r 19 
\ iz 


< 0.6 


45 


RUT nn 




< 47 


< 5.8 




< 52 


< 2.6 


46 


RYLup 




< 37 


< 4.5 




< 23 


< 1.1 


47 


FX T nn 
DA l_,lip 




91 

\ ZJ 


<r- 9 Q 
^ Z.7 




<r- 17 

<-s 1 / 


/ Q 

<v U.7 


48 


Szl02 




<r 17 


<r 9 1 

<^ Z. 1 


970 


27 


1.4 


4Q 


AS 205 




07 


<f 1 1 1 




1 17 

x 11 / 


/ S Q 

<-s J." 


50 


Haro 1-1 




<, o 


^ i i 

V. 1.1 




c -- 7 


^ n i 

\ U.J 


<i 

j i 


Haro 1 -4 




/ Q 


/ 1 1 

<s. 1.1 




^ iq 

\ 17 


/ Q 


S9 


FtnAr 94F 




< 30 


/TO 
< . J.O 




^ R7 


< 4 4 


53 


SR21 




< 120 


< 14.7 




< 19 


< 1.0 


54 


IRS51 




< 32 


< 3.9 




< 9 


< 0.5 


j j 


SR 9 




<r- 91 
s. Zl 


^9 7 
<^ z. / 




^ i/i 
<v it 


/ 7 
<v u. / 


S6 

JO 


V oJ JL/pil 




^ Q 


1 1 

V. 1.1 




c -- q 


< 5 


57 


Rnx4.9r" 




< 13 


< 1.6 




< 7 


< 0.4 


58 


ROX41 A 

1\V J /\ '-r J rt 




< 26 


< 3.2 




< 52 


< 2.5 


SQ 


TR <560 




^ Q 

\ y 


<r 1 9 
<^ 1 .z 




/ -27 


/ 1 c 
<^ 1 .o 


60 


Harn 1-16 
n at u i - 1 u 




<r 18 


< 2.2 




<r 94 

\ ZH- 


<r 1 9 

<-v 1 .Z 


O 1 


Harn 1 1 7 

n ai u i- 1 / 




/ 7 


\ 1 .u 




1 

\ J 


/ o 1 

<v U. 1 


62 


RNO 90 




< 70 


< 8.7 




< 62 


< 3.1 


63 


WaOph6 




< 19 


<2.3 




<25 


< 1.3 


64 


V11210ph 




<33 


<4.0 




<37 


< 1.9 


65 


HD 163296 




< 120 


< 14.4 




<300 


< 15.1 


66 


HD 163296 




< 170 


< 21.1 




<725 


<36.3 


67 


VVSer 




< 150 


< 19.2 




<500 


<25.3 


68 


SSTc2dJ 182900.9 




<53 


<6.5 




<40 


<2.0 


69 


SSTc2dJ182909.8 




<67 


<8.3 




<95 


<4.7 


70 


SSTc2dJ182928.2 




<340 


<42.3 




< 167 


< 8.3 


71 


EC74 




<48 


<5.9 


974 


27 


1.3 


72 


EC82 




< 150 


< 18.1 




< 137 


<6.8 



18 



TAB LE 3 — Continued 





Source 


T b 

^warm 


A/warm^ 


J^waim^ 


7hot b 












(r= 100AU) C 






(/ = 2AU) C 








[K] 


[10 22 cm- 2 ] 


[Mj] 


[K] 


[10 22 crrT 2 ] 


[10" 3 Mj] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


73 


EC90 




< 340 


<42.1 




< 875 


<43.4 


74 


EC92 




< 110 


< 13.7 




< 80 


< 4.0 


75 


CK4 




< 51 


<6.3 




< 32 


< 1.6 


76 


LkFta 348 




<71 


<8.8 




<675 


<33.0 



a 100 K is assumed for the warm component and 1000 K for the hot component. The derived column density and mass of the warm component depend strongly on the assumed 
temperature. A temperature of 150 and 200 K reduces the column density and mass by respectively a factor of ~ 30 and 140. For the hot component 1500 K instead of 1000 K may 
results in a reduction of column density and mass up to a factor of 10. 

b For the unresolved emission a source with r = 50 AU is assumed to obtain an estimate of the column density. 



TABLE 4 

Summary of data for sources with detected H 2 , [Ne II] , or [Fe I] emission 



# Source logL logL logL logL logL T* Age Class log Lx a 

H 2 S(2) H 2 S(3) [Nell] [Fel] IR 
[L ] [L ] [L ] [L ] [L ] [K] [Myr] [ergs" 1 ] 



27 


C7-11 


<-5.6 




<-5.7 


-4.3 


-2.1 


4860 


0.6 


TTs 




21 


Ced 1 10 IRS6 


-5.3 




-5.3 


< -4.8 


-1.8 






1 


28.6<D 


13 


CoKu Tau 3 


<-5.6 




<-5.4 


-4.6 


-2.1 


3800 




CTTs 


30.8< 2 > 


71 


EC74 


-5.0 


-4.0 


-5.3 


< -4.4 


-1.9 






2 


30.5< 3 > 


72 


EC82 


-4.4 


<-3.3 


<-4.4 


< -4.0 


-0.9 


3918 




TTs 




74 


EC92 


-4.8 


< -3.6 


-4.7 


< -4.3 


-1.2 


4273 




TTs 




37 


HD 132947 


<-6.2 


< -5.5 


<-6.4 


-6.3 


-3.2 


10040 




HAeBe 




61 


Haro 1-17 


<-6.0 


<-5.1 


-6.1 


-4.7 


-2.6 


3500 




TTs 




44 


IM Lup 


<-5.4 


< -4.4 


-5.4 


< -5.0 


-1.9 


3918 


0.7 


TTs 




7 


IRAS 03446+3254 


<-5.4 


< -4.4 


-5.2 


< -4.2 


-2.0 






TTs 




17 


IRAS 08267-3336 


<-4.5 


<-3.9 


-4.4 


< -4.0 


-1.1 


5010 


2.0 


TTs 




59 


IRS60 


<-4.9 




-5.2 


< -5.0 


-1.8 










26 


ISO-Cha237 


-5.5 




<-5.2 


< -4.6 


-1.9 


3918 




TTs 


29.2(4) 


3 


Lktfo 271 


<-5.2 


<-4.6 


<-5.2 


-3.8 


-2.3 


4860 




TTs 




55 


SR9 


<-5.2 




-5.5 


< -5.0 


-2.0 


4155 


0.8 


TTs 




70 


SSTc2dJ182928.2 


<-4.2 


<-3.1 


-3.8 




-0.5 










48 


Szl02 


-4.9 


-4.1 


-4.3 


<-4.7 


-2.2 


3918 




TTs 


28.6< 5 > 


42 


Sz73 


<-5.5 


<-4.9 


-5.3 


<-5.2 


-2.1 


3918 


3.7 


TTs 




32 


TCha 


< -6.0 




-5.4 


< -5.6 


-2.8 


5430 


2.1 


TTs 


30.1 <6) 


56 


V8530ph 


<-5.6 




-5.3 


<-5.1 


-2.1 


3800 


0.8 


TTs 




23 


VWCha 


<-4.8 




-4.6 


< -4.5 


-1.6 


4510 


1.2 


CTTs 


31.0< 7 > 


25 


WXCha 


<-5.2 




<-5.3 


-4.2 


-1.9 


4273 


2.3 


TTs 


30.2 (7) 


29 


XX Cha 


<-5.7 




-5.5 


-4.5 


-2.2 


3800 


14.6 


TTs 





a X-ray luminosities taken from: (1) Carkner et al. (1998), (2) Audard et al. (2006), (3) Preibisch (2003), (4) Feigelson & Lawson (2004), (5) Gondoin (2006), (6) Alcala et al 
(1997), (7) Stelzer et al. (2004). 
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